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ABSTRACT 

If 7-ray bursts arc cosmological in origin, the sources are expected to 
trace the large-scale structure of luminous matter in the universe. I use a 
new likelihood method that compares the counts-in-cells distribution of 7-ray 
bursts in the BATSE 3B catalog with that expected from the known large-scale 
structure of the universe, in order to place a constraint on the distance scale 
to cosmological bursts. I find, at the 95% confidence level, that the comoving 
distance to the "edge" of the burst distribution is greater than 630 h"^ Mpc 
{z > 0.25), and that the nearest burst is farther than 40 Mpc. The median 
distance to the nearest burst is 170 Mpc, implying that the total energy 
released in 7-rays during a burst event is of order 3 x 10^^ h^"^ ergs. None of 
the bursts that have been observed by BATSE are in nearby galaxies, nor is a 
signature from the Coma cluster or the "Great Wall" likely to be seen in the 
data at present. 

Subject headings: gamma rays: bursts — large-scale structure of universe — 
methods: statistical 
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1. Introduction 

The origin of 7-ray bursts is still unknown and is currently the subject of a "great 
debate" in the astronomical community. Do the bursts have a Galactic origin ([Lamb 1995|) 
or are they cosmological ([Paczyhski 1995| )? And what is their distance scale? 



In this Letter, I do not attempt to answer the first question, but rather, I show that 
if one assumes that 7-ray bursts are cosmological in origin, one can begin to answer the 
second question and place a constraint on the distance scale to the bursts. This is because 
cosmological bursts are expected to trace the large-scale structure of luminous matter in the 
universe (|Lamb fc Quashnock 1993| , hereafter, LQ). The constraint comes from comparing 



the expected clustering pattern of bursts on the sky — which will depend on their distance 
scale because of projection effects — with that actually observed. The observed angular 
distribution is in fact quite isotropic ( Priggs et al. 199^ ); hence, only a lower limit to the 
distance scale can be placed because a sufficiently large distance will always lead to a 
sufficiently isotropic distribution on the sky. 

Hartmann and Blumenthal (1989) first used the absence of a significant angular 
correlation function in a small sample of 7-ray bursts to put a lower limit of 71 Mpc 
on the distance scale to the bursts.^ This lower limit was subsequently improved to 
150 h-^ Mpc ([Blumenthal, Hart mann, fc hinder 1994[) using the larger BATSE IB catalog 
( [Fishman et al. 1994[) . 



Here I use a powerful new likelihood method, which I had previously developed to 
analyze repeating of 7-ray bursts in the BATSE IB and 2B catalogs ( [Quashnock 1995| ), to 
compare the observed counts-in-cells distribution in the new BATSE 3B catalog ([Meegan[ 



[ % follow the usual convention and take h to be the Hubble constant in units of 100 
km s~^ Mpc~^. 
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et al. 1996) with that expected for bursts at cosmological distances. I describe this method 
and calculate the expected counts-in-cells distribution in § 2.1, using the angular correlation 
function computed in § 2.2. I present my results on the burst distance scale in § 3, and 
discuss some implications of this work in § 4. This work was presented in preliminary form 
elsewhere ( IQuashnock 1996| ). 

2. Likelihood Method 

2.1. Counts— in— Cells Distribution 

Let A^^ceii be a large number of circular cells, each centered on a random position on 
the sky. Each cell is of fixed solid-angle size f2 = 27r(l — cos^rad), where 6'i.ad is the angular 
radius of the cell. I set the number of cells to be such that any part of the sky is covered, 
on average, by one cell; hence, Nceii = 47r/f2. Let Cat to be the number of these cells 
having 7-ray bursts in them, out of the A^tot = 1122 in the BATSE 3B catalog, where 
A^ = 0, 1, 2, . . .. Then I define the observed counts-in-cells distribution, 

= CM/N,,n , (1) 

which is the probability that a randomly chosen cell of size Q has A^ bursts in it. The 
counts-in-cells distribution contains information about clustering of 7-ray bursts on scales 
comparable to the angular size 6'rad of the cell. 

I now define Qat to be the counts-in-cells distribution that is expected if 7-ray bursts 
are cosmological in origin and trace the large-scale structure of luminous matter in the 
universe. This expected distribution depends on only one unknown parameter, the effective 
distance D to 7-ray bursts (which I define below), because the angular clustering pattern 
of bursts on the sky will depend by projection on this distance. 
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The likelihood C measures how hkely it is that the observed count s-in-cells distribution 
P/v is drawn from the expected distribution Q^. Since Qn depends on the unknown 
effective distance D to 7-ray bursts, the hkehhood is really a measure of how likely a given 
value of D is. It is given by the multinomial expression 

^ = Tr^,IlQN^- (2) 

Combining equations (1) and (2), and noting that the factorial term is parameter 
independent, I obtain 

log £ = A^ccu ^PN^ogQN + constant . (3) 

The expected counts-in-cells distribution Qn is most easily computed using the 
generating function Q(A) = J2nQn^'^ ( Peebles 1980| ; Palian fc Schaeffer 1989|) and writing 



it in terms of the means, over the cell size, of the irreducible A^-point angular correlation 



functions wn ([White 1979| ; [Balian fc Schaeffer 1989| ): 

/N\^ ..1 

(4) 



Q(A) = exp 



00 /ivr\N 
Af=2 



where the density (A^) = A^tot/^ceii is the mean number of bursts in a cell. Indeed, the 
Poisson distribution, Qn = (A^)^exp {—{N)) /N\, is obtained (by comparing powers of A) 
when all correlations Wn are zero (uniform distribution) and only the first term in the sum 
in equation (4) is included. 

Now for cosmological 7-ray bursts, the clustering is expected to be weak (LQ) and can 
be adequately described by one number, W2, the mean of the two-point angular correlation 
function w{6) over the cell: 

W2 = ^ f dn^dn2wi9i2) . (5) 



fi2 

Then only the first two terms in the sum in equation (4) contribute to the generating function 
and the counts-in-cells distribution is a convolution of two terms: Qn = Y.m=o ^n-m Bm, 
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where the first is a Poisson uncorrelated background of single bursts of density (A^) — {N)'^W2, 
and the second a uniform smattering of correlated pairs, with density {N)'^W2/2, in which 

(6) 

Once the angular correlation function is known, W2 can be computed from equation 
(5), and the expected count s-in-cells distribution Qm can be found from equation (6). The 
likelihood C can then be computed from equation (3), and it is a function of the unknown 
effective distance D through the dependence of W2 on D. 



2.2. Angular Correlation Function 

For simplicity, I assume that r2o = 1 and A = 0, and that the large-scale structure 
clustering pattern is constant in comoving coordinates. The results (cf. § 3) are, in fact, 
insensitive to these assumptions because of the small redshifts that are involved. 

The Limber equation ( Peebles 1980| ), which relates the angular correlation function 
w{9) to the spatial one, can be recast (Peacock 1991| ; LQ) in terms of the dimensionless 
power spectrum A^(/c): 

^ ir dr 7rA^(fc) Jo(fcrg) dk/k' 

^ [/crr20(r)rfrf ' 

Here the selection function 0(r) is the probability that a source at comoving distance 
r produces a burst that is in the BATSE catalog, i.e., a burst with apparent fiux greater 
than the limiting fiux of the BATSE survey. The selection function is a convolution of a 
source function, representing the luminosity function and the spatial density of the sources. 
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with an observer function, representing the efficiency with which BATSE detects bursts of 
a given flux. 

Now the cumulative Cmax/C*min distribution of 7-ray bursts seen by BATSE begins 



to roll over from a —3/2 power law for bursts fainter than Cmax/C'min ~ 10 ( |Meegan~et 



al. 1992|) . Since this is many times above threshold, it suggests that BATSE sees most of 



the source distribution and that this distribution is not spatially homogeneous. To the 
extent that this is the case, the observer function is unity for nearly all the burst sources, 
irrespective of brightness, so that BATSE sees nearly all of the burst sources. The selection 
function (f){r), then, depends only on the luminosity function and the mean space density of 
the sources. 

I define D as the distance beyond which 0(r) drops appreciably; thus, D is the 
effective distance to the "edge" of the source distribution in the BATSE catalog. One could 
approximate the selection function to be unity out to a comoving distance and zero 
beyond, but evolution of the luminosity function and/or the mean space density of the 
sources will change the form of 0(r) from 9{r — D); more generally, I define the effective 
distance D as ^D^ = J (f){r) dr (LQ). Thus, D is not the distance to the very dimmest 
burst in the BATSE catalog, but rather the typical distance to most of the dim bursts in 
the sample. 



Equation (7) then becomes, upon rescaling y = r/D, 

HO) = — yUy A\k) UkyD9) dk/k' 
D Jo Jo 



so that the angular correlation function satisfies the well-known scaling relation ([Peebles 
T980D w{9) = D-^W{D9) as a function of D. 



I take the power spectrum that characterizes the large-scale clustering of 7-ray burst 
sources to be the same as that determined from a redshift survey of radio galaxies ( Peacock 
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and Nicholson 1991): 



A2 



(fc) 



0.129 (k/k^Y 



(9) 



where = 0.025 /i Mpc'^ This power spectrum is characteristic of moderately rich 
environments, and is intermediate between that of ordinary galaxies and clusters (LQ). 
Because the exact bias factor relating the clustering of 7-ray burst sources to that of 
luminous matter is unknown, such an intermediate Ansatz is reasonable. In any case, the 
resultant distance limit depends only weakly on the bias factor (roughly as the square root). 

Substituting equation (9) into equation (8), one obtains the expected intrinsic angular 
correlation function w{6) of 7-ray bursts as a function of their effective depth D. Figure 1 
shows w{e) {dashed lines) for D of 500, 1000, and 2000 h'^ Mpc. Note the scaling of w{e) 
as a function of D. 

Now the observed angular correlation function, wlO), is smeared at small angular 
scales because of finite positional errors. Each burst in the BATSE catalog is assigned a 
positional uncertainty 6'err corresponding to a 68% confidence that the true burst position is 
within an angle 6'err to the position listed in the catalog. Defining = (1 — cos(6'err)) /1-14, 
the smearing function can be approximated as a Gaussian of the angular separation 
between the true and observed positions of the burst: dP/dQ = exp(— A0^/2(T^)/27rcr^. If 
all burst positions are smeared by the same amount, cr^, the observed angular correlation 
function is a convolution of the intrinsic correlation function with the smearing function 
and the modified Bessel function Jo ( Partmann, hinder, fc Blumenthal 1991| ): 



Figure 1 shows w{6) {solid lines) for D of 500, 1000, and 2000 Mpc, with smearing 
of burst positions of ^crr = 3?8, the median value in the BATSE 3B catalog. Note the 
smearing of the correlation function on scales smaller than 6'err, along with the aliasing 




(10) 
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of some small-scale power to angular scales comparable to the positional smearing. It is 
this smeared correlation function that I use in equation (5) when calculating W2 and the 
expected counts-in-cells distribution Qn/vh calculating the likelihood £ as a function of D. 

The cell size ^rad is chosen in order to maximize the sensitivity of detection, or signal- 
to-noise ratio, given the strength of the signal expected. The signal S is given by the total 
number of correlated pairs in the cells (cf. above eq. [6]): S = Nceii{N)'^W2/2. The noise N 
is the square root of the total number of pairs ( Peebles 1980| ; LQ): = {Ncc\\{N)'^ /2)^^'^ . 
For a sample of 1122 bursts (corresponding to the total number of bursts in the BATSE 3B 
catalog) with positional smearing of 6'crr = 3?8, the signal-to-noise ratio is maximized when 
cells of ^^rad = 5° are used. 

3. Results 

Figure 2 shows the likelihood of the BATSE 3B catalog data as a function of the 
effective comoving distance D, calculated using cells of size ^rad = 5°. The likelihood is 
normalized to that expected for an isotropic distribution on the sky. At large values of D 
(the maximum value allowed is D = = 6000 /i^^ Mpc, the size of the horizon in a closed 
universe), the likelihood goes to unity, because by projection a sufficiently large distance 
will always lead to an isotropic distribution on the sky. Note also that there is no value of 
D for which the likelihood is greater than 1; thus, the maximum likelihood value for D is 
i?H, and the 3B data are consistent with isotropy. 

The solid line in Figure 2 shows the likelihood for a positional smearing of ^err = 3?8, 
corresponding to the median value in the SB catalog. To illustrate the dependence of these 
results on positional errors, I also show {dashed line) the results for a larger positional 
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smearing of ^err = 6?6 (with cells of size ^rad = 9° to maximize the signal-to-noise ratio). 

Small values of the effective comoving distance to 7-ray bursts are unlikely, according 
to Figure 2; I find, at the 95% confidence level, that for the 3B median positional error of 
6'err = 3?8, D must be greater than 630 h'^ Mpc, corresponding to a redshift z > 0.25. If 
the positional errors are larger than quoted and are better characterized by Ocrr = 6?6, these 
results are only slightly weakened; at the 95% confidence level, D must be greater than 
500 h^^ Mpc, corresponding to a redshift z > 0.19. 

These limits are not sensitive to earlier assumptions (§ 2.2) on cosmology and clustering 
evolution, since these only become important at higher redshifts. They are also conservative 
limits, in that a constant median value for the positional errors was used rather than the 
entire distribution of errors. This is because the bright bursts, which ostensibly are nearer 
to us, are more clustered (by the scaling property of the correlation function, eq. [8]) and 
are responsible for the bulk of the expected signal, but, in fact, have smaller errors than 
the median value. The faint bursts, which are far away, are hardly clustered to begin with 
(even before smearing), but have errors larger than the median value. Hence, the expected 
clustering pattern has been smeared more by using a constant median value (this permits 
a simpler calculation and allows eq. [10] to be used) than by smearing using the entire 
distribution of errors. Therefore, the counts-in-cells statistic has been weakened somewhat, 
and thus the quoted lower limits are, in fact, conservative. 



■^Graziani & Lamb (1996) compare the 3B positions with those from the IPN network and 
conclude that the systematic errors are larger than the 1?6 value quoted in the 3B catalog. 
Their best-fit model gives a median positional error of 6?6. 
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4. Discussion 

If 7-ray bursts are cosmological and trace the large-scale structure of luminous matter 
in the universe, and their positional errors are as quoted in the 3B catalog, then the lack 
of any angular clustering in the data implies that the observed distance to the "edge" of 
the burst distribution must be farther than 630 Mpc. Since there are 1122 bursts in 
the catalog, an effective limit on the nearest burst to us can be placed by convoluting the 
likelihood as a function of D (Fig. 2) with the nearest neighbor distribution of 1122 bursts 
inside a sphere of radius D. I find that the nearest burst must be farther than 40 Mpc 
at the 95% confidence level, and farther than 10 h^^ Mpc at the 99.9% level. At this level 
of confidence, then, none of the bursts that have been observed by BATSE are in nearby 
galaxies. A signature from the Coma cluster or the "Great Wall" (~ 70 Mpc) is not 
likely to be seen in the data at present, since only a few bursts could have originated from 
these distances. Indeed, a search for such a signature ( [Hartmann, Briggs, fc Mannheini 



1996|) found no compelling evidence for anisotropy in supergalactic coordinates. 



The median distance to the nearest burst is 170 Mpc. Since the brightest burst in 
the 3B catalog has a fluence of 7.8 x 10~^ ergs cm~^ in 7-rays, this implies that the total 
energy released in 7-rays during a burst event is of order 3 x 10^^ ergs. 

As the number of observed 7-ray bursts keeps increasing, the distance limit will 
improve. In fact, LQ showed that, with 3000 burst locations, the clustering of bursts might 
just be detectable and would provide compelling evidence for a cosmological origin. If it is 
not detected, the redshift to the "edge" of the bursts would be put at z ~ 1 or beyond. 



I would like to acknowledge useful discussions with Carlo Graziani, Don Lamb, Cole 
Miller, and Bob Nichol. This research was supported in part by NASA through the Compton 
Fellowship Program — grant NAG 5-2660, grant NAG 5-2868, and contract NASW-4690. 
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6 (degrees) 

Fig. 1. — Angular correlation function and its dependence on the effective comoving distance 
D to 7-ray bursts. Shown are both the intrinsic, w{6) {dashed line), and smeared, w{9) {solid 
line), correlation functions, in decreasing amplitude, for D = 500, 1000, and 2000 h^^ Mpc. 
The smearing corresponds to — 3?8, the median value in the BATSE 3B catalog (see 
text) . 
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Fig. 2. — Likelihood of the BATSE 3B catalog data as a function of the effective comoving 
distance D to 7-ray bursts, shown with a smearing of ^err = 3?8 and 6?6 (see text). 



